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Abstract 


Introduction. Slender structures of subsea energy production systems are under constant influence of currents and waves. 
Hydrodynamic loads result from the interaction of subsea pipelines, umbilicals, equipment supports with fluid flows, and 
lead to the vortex formation in the area behind the structures. Vortex-induced forces are the sources of the cyclic loading. 
They accelerate gradually the fatigue damage, which may result in a failure. One of the ways to reduce the loads on subsea 
structures is to alter the shape of a cross-section, taking into account the flow regime. Dependence of the resulting 
hydrodynamic loads on the cross-sectional shape and relative position of structures has not been studied in details for the 
uniform flow in the critical mode. The current work is aimed at filling this gap. The research objective is to consider the 
impact of the distance between the structures, and also, the presence of a D-shaped structure, placed upstream relative to 
the group of three cylinders of different cross-sectional shapes. 

Materials and Methods. The computational fluid dynamics approach was used in this work for numerical simulations of 
vortex-induced forces in the ANSYS Fluent software for cylinder with D = 0.3 m. Modelling was conducted with the 
Detached Eddy Simulation (DES) method, which combined advantages of the Reynolds-averaged Navier-Stokes equation 
(RANS) method and the Large Eddy Simulation (LES) method. The object of the research was the system of four 
structures in the 2D computational domain, which included the upstream D-shaped cylinder and the main group of three 
cylinders with the circular, squared and diamond shapes of the cross-section. The transient process was considered, where 
structures were under the influence of the uniform flow in the critical regime at Re = 2.5x10°. 

Results. Five sets of data were obtained in simulations for the time-dependent coefficients of the lift and drag forces: for 
the main system — of the D-shaped, circular, square and diamond structures, and also for the four systems — of only 
D-shaped, only circular, only square and only diamond shaped structures. Additional analysis was conducted for the effect 
of the distance between the structures on the amplitude of fluctuating hydrodynamic force coefficients. The obtained 
results are presented as time histories of coefficients of the lift and drag forces, frequency analysis and contours of 
velocity, pressure and vorticity fields. The results indicate a positive effect of the upstream D-shaped structure on reducing 
the drag force, acting on the central structure in the group of three cylinders located downstream. 

Discussion and Conclusion. The results of the performed studies facilitate the informed decisions regarding the 
arrangement of subsea structures in a group of four objects, depending on the cross-sectional shape and the distance 
between the structures. The upstream D-shaped structure provides reducing the hydrodynamic drag force acting on the 
central structure in the downstream group of three structures, thereby slowing the fatigue accumulation and increasing 
the time of safe operation. 
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AHHOTalna 

Beedenue. J\nmuusie u y3Kve B MOMepeyHHKe KOHCTpyKUMM MOPCKHX IHEprosOOKIBalOWHX CHCTeEM HaxXOATCA MO T0- 
CTOAHHBIM BO3eMCTBHeM TeyeHH wu BOTH. [uapodqMHaMMyecKne Harpy3KH ABJIAIOTCA Pe3yJIbTATOM B3aMMOJeHCTBUA 
MOABOAHBIX TpyOOMpoBOOB, WaHToKaOeseH, oop OOOpyAOBaHHA C MOTOKOM 2KMAKOCTH HM MPHBOAAT K OOpa3z0BaHHto 
BUXpel B 30He 3a KOHCTpyKuMaAMu. BuxpeoOpa30BaTeJIbHble CHJIbI CILYKAT HCTOUHHKOM UMKIMYeCKOrO Harpy2KeHHA UH 
TIOCTEMCHHO YCKOpAIOT YCTaIOCTHOe pa3spylIeHHe, YTO MOXeT MpHBeCTH K aBapHAM. OJHUM U3 CHOCOOOB CHH>KeHHA 
Harpy30K Ha NOABOJHbIC KOHCTPYKUMH ABJIACTCA H3MCHEHHe (POPMBI HX MOMepe4yHoro CeyeHHA C YUCTOM pexKHMa HOTOKAa. 
HeyoctaTouHo u3y4eHo, KaKHM 06pa30M HTOTOBbIe rH pOMHAaMMYeCKHe Harpy3KU 3aBHCAT OT (POPMBI WoMepeyHoro ce- 
YeHHA HW B3AMMHOTO paciOJIOXKeCHHA HA3BAHHBIX BbILIe 3ICEMCHTOB CHCTeM, HaXOJ{ALIMXCA B paBHOMEPHOM KPHTH4YeCKOM 
mioToke. IIpeyctapyeHHad Hay4ynad padoTa MpH3BaHa BOCHOJHUTS ITOT Mpoben. Lemb ucceqyoBaHHA — PpaccMOTPeTb B 
JJaHHOM KOHTeKCTe 3Ha4eHHe paccTOAHHA MeX*KAy KOHCTpyKUMAMH, a TakKxKe HaMuve NosyKpyrion D-oOpa3sHoi KOH- 
CTPpyKIUMH, pa3sMeléHHOH Nepey rpynioli 43 Tpéex WWM APOB C pasHbIMM MonepeuHbIMM CeyeHHAMH. 

Mamepuanoi u memooot. Jna wicneHHoro MoyeuMpoBaHHA BUXpeoOpa3z0BaTeJIbHbIX CHI HCIHOIb3OBaJICA MCTO], BbIMHC- 
JMTeIbHOM JHHAaMUKH (pIOMAOB B IporpamMMe ANSYS Fluent aa uA ApoB AMametpom D = 0,3 m. MoyemupoBanue 
BBIMOJIHCHO Me€TOJOM HelIpHcoeqHHéHHBIX BUxpeli DES, koTopblii coueTaeT B cebe NpeMMylecTBa MeTOa ycpeqHeH- 
Horo 10 PeitHombacy ypaBHenua Happe-Ctoxca RANS nu metoyia kpynHbIx Buxpeli LES. B kayectBe oOnexta uccieqo- 
BaHHA pacCMaTpHBallach CHCTeMa, COCTOAMAd W3 YeTLIPEX KOHCTPYKIMM B BbIYMCIMTeIbHOM JOMeHe B 2D, BKuIOUAA 
CTOALIMM BbIMe 10 TEYeEHHIO NOMYKpyribli WWIMHAp HU OCHOBHY!O rpyninly 43 Tpéx WHIMHApOB Kpyrioi, KBaypaTHoH u 
POMOoOBH HOM (POpMbI WoMepeyHoro CceyeHHA. ITH KOHCTPYKUMU B yCIOBHAX H€YCTAaHOBHBIIeroca polecca HaxOWATCA 
lO AelicTBHeM paBHOMepHOro NoTOKa KpuTMYecKOro pexkuma pH Re = 2,5x10°. 

Pe3yibmamot ucciedosanua. B pe3yibTaTe MOJeIMpOBaHHA MOJYYCHI MATb HAOOPOB AaHHBbIX JIA M3MCHAIOLIMXCA BO 
BpeMeHH KOOPPUUMEHTOB BUXpeoOpa3z0BaTeNbHBIX MOAbEMHOM CHJIbI HM CHJIbI COMPOTHBIICHHA: JIA OCHOBHOM CHCTeMBI 
M3 MOYKpyrsiol, Kpyriol, KBadpaTHOM H pOMOOBHAHOM KOHCTpyKUMH, a TakoKe JIA YeTbIPEX CHCTEM 3 TOJIbKO TIOJy- 
KPYTJIBIX, TOJIBKO KPYTJIBIX, TOJIBKO KBaPaTHBIX H TOJIbKO pOMOOBHHBIX KOHCTpyKUHH. JJOMOJHUTeIbHO IpOBeyGH aHa- 
JIM3 BIMAHHA paCCTOAHHA MOK Ly KOHCTPyYKUMAMH Ha aMIWVIMTyAy KoeOaHnit KoIpPUUHeHTOB rH pOAMHAaMHYeCKHKX CHIL. 
Tlomy4enHble pe3yIbTaTbl MpeACTaBIeHbI B BUC KOIPPUUMEHTOB NOAEMHOM CHJIbI M CHJIbI COMpPOTHBIICHHA B MHa- 
MHKe, aHasI43a YaCTOT HM KOHTYPOB MOJIeH CKOPOCTH, WaBIeHHA, 3aBUXpEHHOCTH. Pe3ybTAaTbI NO3BOJIAKOT YCTAHOBUTb 
HOJIOXKUTEIbHOE BIMAHME CTOALMICH BbIMEe MO TEYCHHIO NOUYKpyrsiol KOHCTpyKWHH Ha CHW KeHHE CHJIbI CONPOTHBIICHHA 
Ha WeHTpasIbHy1O KOHCTPyKIHH0 B rpylitie 43 TpéxX WHIMHApOB HDKe 110 TeYeHHIO. 

O6cystcoenue u 3akitovenue, Pe3yibTaTbI IpOBeCAEHHBIX HUCCIeTOBaHHM MO3BONAIOT IPHHAMaTb OOOCHOBAHHEIEe pellie- 
HUA [IX PaCCTAaHOBKM MOPCKHX KOHCTpyKUMM B pyle 43 YeTLIPEX OOLEKTOB B 3ABHCHMOCTH OT (POpMBI HoMepeyHoro 
ce4ueHHaA HU paccTOAHMA MexKy HMMM. YcTaHOBKa NOYKpyrioH KOHCTPyKUMH BIE TO TEYeHHIO NO3BOJIAeT CHH3HTB 
TH pOAMHAMHYeCKy!O CHJIY COMPOTHBJICHHA Ha WeHTpasIbHy!O KOHCTPyKUMHO B rpyniie 43 TpéxX KOHCTpYKUMH HWKe 110 
TeC4eHHIO, YTO 3aMeIAeT C€ YCTANOCTHOE paspylleHHe HW YBeJIMYMBAeT CPOK SKCIIyaTauHu. 


Ko1r04ueBble CJIOBa: BHxpeoOpa30BaTelIbHble CHIBI, Kospunnent JIOOOBOTO CONpOTHBJICHHA, Ko3dulnent T1IObEMHOI 
CHIBI, paBHOMepHBIit TLOTOK 
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Introduction. Operation and construction of modern offshore systems, specializing on the energy production, 
extraction of resources or carbon capture and storage, require evaluation of the impact of environmental flows on 
equipment and structures. An increased fatigue failure in subsea structures, such as pipelines, risers, cables, piles, 
equipment supports, may come from the vortex shedding phenomenon. The problem is particularly important when 
slender structures are designed to reach deep waters to connect subsystems together. The layout of subsea systems comes 
with the arrangement of structures with different geometry, hydrodynamic properties, and their position in proximity to 
each other. The interference of wakes from these structures and vortex formation patterns is sometimes challenging to 
predict at very high Reynolds numbers due to the turbulent nature of the flow. 

Differences in the flow over a standalone cylinder and two cylinders in tandem are discussed in [1], and three vortex 
shedding regimes are identified for tandem structures. These vortex formation modes include the extended-body regime 
at 1.0 < L/D < 1.8, where L/D is spacing ratio, commonly used to quantify the distance between centres of neighbour 
structures. So that, Z corresponds to the distance, and D is the diameter of the structure. In [1], increase of the spacing 
ratio to 1.8 < L/D < 3.8 leads to the reattachment regime, where shear layers detach from the upstream structure and 
reattach to the front side of the downstream structure, so that vortices are formed behind this downstream object. Further 
growth of the spacing ratio, above L/D > 3.8, introduces the co-shedding regime, where a separate vortex is formed from 
the upstream structure and from the downstream structure. Another fundamental research investigated vortex dynamics 
in details through experimental research [2]. One of the following fundamental studies [3] experimentally investigated 
the vortex shedding frequencies of two staggered identical circular cylinders with the Reynolds number Re varying from 
3.2x10* to 7.4x104, and two fixed side-by-side cylinders at the Reynolds number of 2.5x10* were earlier considered in [4]. 
These investigations provide an important foundation for modern studies in terms of the known effects in fluid forcing 
and vortex shedding patterns. 

Significantly more recent investigations are performed numerically [5] to study the effect of spacing on loads and 
vibrations for two tandem cylinders at subcritical Reynolds numbers, and for specific cases, like a group of mixed large 
and smaller structures [6]. The latter work [6] numerically investigates fluid force coefficients and observes the vortex 
formation pattern on three identical rigid circular cylinders in proximity to a square cylinder. A parametric study is 
conducted in [7] for three identical stationary circular and D-shaped cylinders placed close to a square cylinder at 
Reynolds number 3900 in both linearly and parabolic sheared flows. 

Considering the impact of cross-sectional shapes further, a numerical study is conducted by [8] for a flow over six 
identical stationary cylinders having different cross-sectional shapes at Reynolds number of 2.5x10° in the uniform and 
linearly sheared flow. Rectangular cylinders are investigated in details in [9, 10, 11], where one of the most impactful 
factors for hydrodynamic loads is the aspect ratio of rectangle sides. The works [9, 10] provide new experimental data 
and attempt to develop semi-empirical methods of predicting the response of structures. Further steps in improving the 
modelling approaches for the structural vibration of rectangle-shaped objects under the hydrodynamic excitation are 
performed in [11]. Another branch of studies considers a flow over a sphere [12, 13, 14], while still leaving issues of the 
impact of the cross-sectional shape of subsea structures open. 

Further research on diverse cross-sectional shapes is performed in [15, 16, 17], where triangular and diamond cross- 
sections are studied in comparison. Research [15] is focused on the sensitivity to the corner sharpness for the diamond 
(rhomb) cross-section. Work [16] investigates effects from diverse cross-sections for the system with a rotational degree 
of freedom, when subjected to flow-induced vibration, and study [17] uses cross-sectional shapes for the energy 
harvesting with fluid-structure interaction. 

Following published results, the current work aims to investigate the drag reduction when three structures of different 
cross-sectional shapes are located around a circular cylinder to observe the wake interference and the vortex formation 
pattern between these structures for the spacing ratio L/D varying from 1.67 to 2.83 in the uniform flow at the Reynolds 
number of 2.5x10° with the computational fluid dynamics approach. Specifically, the drag reducing ability of the upstream 
structure is of the research interest. The considered layout is a combination of cylinders placed in tandem, side-by-side 
and staggered position, with a mixture of cross-sectional shapes. 

Materials and Methods. A system of four cylinders with a diameter (characteristic size) of 0.3 m with different cross- 
sectional types is considered in this study. The selected cross-sectional shapes and positions of structures are shown in 
Figure 1. Fluid forces and the flow interference are studied for the spacing ratio of L/D varying from 1.67 to 2.83 in the 
uniform current at the Reynolds number of 2.5x10° and the corresponding velocity of 0.837 m/s. The study focuses on 
the specific layout, where structures are positioned relative to each other in a mixed tandem (cylinders | and 3), side-by- 
side (cylinders 2 and 3, and cylinders 3 and 4) and staggered (cylinders | and 2, and cylinders | and 4) configuration. 
Cross-sectional shapes considered include half-circle, square, circle, and diamond. 
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Periodic 


Shadow 


Fig. 1. Considered structures in the computational fluid domain 


Computational fluid dynamic simulations are performed in the domain with a size of 47D x 27D. The top and bottom 
boundaries are located at a distance of 13.5D away from the center of the circular cylinder, periodic and shadow properties 
are assigned to these boundaries. The left boundary serves as a velocity-inlet, which is located at the distance of 20D from 
the centre of the circular structure 3 in the domain. The value of gauge pressure is set to zero at the pressure-outlet set at 
the right boundary. No-slip conditions are applied to the cylinders. 

The flow around cylinders is simulated using the computational fluid dynamics software ANSYS Fluent, where the 
finite volume method is implemented to solve the Navier-Stokes system. The incompressible flow is considered, and the 
2D DES transient simulations are conducted with the A-w SST turbulence model. Time integration is performed using 
the second-order implicit transient formulation with a time step of 0.01 s, and the PISO algorithm is used as the solver. 

The DES approach connects capabilities of the Reynolds-Averaged Navier-Stokes (RANS) and Large Eddy 
Simulation (LES) methods [8]. The RANS governing equations for the incompressible flow are as follows: 


8(pu; 
= =0, (1) 
A(pu;) 9 ,— —— _ op , Oy p 
A #3 (Puie tenia = <P Ox) (2) 


where p is mean pressure, u, is average Cartesian components of the velocity vector, pu;u; are Reynolds stresses, p is 


density of the fluid, and t,; is mean viscous stress vector components, which could be expressed as: 


= OU 
Oty = {% al 3) 


where pl is dynamic viscosity. 
The Large Eddy Simulation (LES) system of equations for the incompressible flow can be written in the following way: 


du; 

— =0, 4 

a, (4) 
i, D 24 Oty 
sca )= 1 op, Ou; un (5) 
ot Ox, pox; Ox,Ox; Ox, 


where U; and P represent the resolved filtered velocity and pressure, respectively. 
The diffusion term of the DES model is given by 
Y, =pB'koF pes, (6) 
where " is a constant, k stands for fluctuation of the turbulent kinetic energy, w is specific energy dissipation rate, and 
F pgs is as follows: 


Foes = mas a, ‘} (7) 


des max 
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where Cas is a constant, Ajax is local maximum grid map A = (Aj, Ao, A3)!*. Further, L; is turbulent length scale: 


Je 


= oe 8 
(Bo (8) 
The DES-SST model uses the following zonal formulation: 
L 
Foes = mas| : (1 Fsr).t} (9) 
des‘ max 


where F'ssr = 0, F1, Fo, and F), F are mixed functions of the SST model. 

Table 1 below provides results of the mesh independence test for the uniform flow of Re = 2.5105. The mesh settings 
are adopted from [7], and the accuracy of the grid is demonstrated by comparisons in Table 1. All subsequent analysis in 
this paper is performed with Mesh 2, and the results include signals and frequencies of the fluid force coefficients and an 
indication of the vortex shedding pattern. The drag force fluctuations are presented in this work in terms of the drag force 
coefficient Cp which comprises the mean drag coefficient Cpo and the fluctuating drag coefficient Cp", as follows: 


Che Cages. (10) 


The lift force fluctuations are presented using the lift coefficient C;. 


Table 1 
Mesh independence test results 
Cases Cbo Number of cells Strouhal number 
Re = 2.5x10° 
Current study 
Mesh 1 0.98 63,345 0.24 
Mesh 2 1.08 86,478 0.24 
Mesh 3 1.08 131,041 0.24 
Published data 
Lehmkuhl, et al. (2014) (LES) [18] 0.833 - 0.238 
Achenbach&Heinecke (1981) (Experiment) [19] 1.135 - 0.230 
Re = 3,900 
Current study, Mesh 2 0.93 86,478 0.18 
Wornom, et al. (2011) (VMS-LES) [20] 0.99 = - 
Re = 3.6x10° 
Current study, Mesh 2 0.4100 86,478 = 
Porteous, et al. (2015) (URANS) [21] 0.4206 - =: 
Nazvanova, et al. (2022) (URANS) [22] 0.4657 74,496 = 


Results. In this study, simulations are performed in two series. The first series is focused on recognising the overall 
effect of various cross-sectional shapes, placed at L/D = 2.00 from each other. Cylinder numbers here correspond to the 
ones used in Figure 1, according to the cylinders position. The calculation results for this set are reported in Table 2, in 
comparison to the case of structures with mixed cross-sections. The second series provides an insight into the impact of 
L/D ratio on hydrodynamic loads observed for the mixed cross-section case only, as in Figure 1. These results are 
summarised in Table 3 and allow defining the drag reducing effect of the D-shape upstream structure on loads when 
cross-sections are different. 
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Simulation results for the same arrangement with different cross-sectional shapes for L/D = 2.00 


Table 2 


LID =2.0 ei oi nie ia All circular | All square All D-shaped All diamond-shaped 
Groseeectons in Riga structures structures structures structures 
Cylinder 1 
Coo 0.45 0.28 0.49 0.48 0.46 
ch 0.19 0.14 0.43 0.28 0.17 
CL 0.18 0.04 0.07 0.13 0.28 
Cylinder 2 
Coo 0.87 0.39 0.89 0.64 0.70 
cf 1.03 0.38 1.38 0.49 0.60 
Cr 1.38 0.53 0.07 0.11 0.07 
Cylinder 3 
Co 0.3 0.24 0.26 0.35 0.70 
cE 0.45 0.37 0.91 0.68 0.63 
CL 1.03 0.73 1.26 0.30 0.73 
Cylinder 4 
Cpo 0.98 0.39 0.83 0.61 0.99 
cf 0.90 0.27 0.96 0.44 1.05 
Cr 1.22 0.81 1.63 0.37 1.11 
Table 3 
Simulation results for the mixed cross-sections at various L/D ratio 
Cylinder 1 Cylinder 2 Cylinder 3 Cylinder 4 
L/D 
Co | CH CL Co | Ch CL Co | Ch Cr Co | Ch Cr 
1.67 0.14 0.18 0.15 0.98 1.17 1.44 0.48 0.57 0.66 0.93 0.75 1.13 
1.83 0.33 0.21 0.18 0.93 1.54 1.40 0.33 0.50 0.86 1.01 1.36 1.15 
2.00 0.45 0.19 0.18 0.87 1.03 1.38 0.30 0.45 1.03 0.98 0.90 1.22 
2.17 0.49 0.05 0.19 0.81 0.14 1.11 0.25 0.06 0.83 0.96 0.38 1.10 
2.33 0.37 0.12 0.00 0.78 0.34 1.03 0.25 0.06 0.28 1.00 0.07 0.79 
2.50 0.57 0.16 0.20 0.76 0.89 1.18 0.26 0.47 0.77 1.01 0.64 1.06 
2.67 0.57 0.02 0.27 0.78 0.16 1.37 0.28 0.08 0.82 1.02 0.03 0.96 
2.83 0.58 0.21 0.21 0.84 0.91 1.70 0.20 0.52 1.16 0.93 0.68 1.03 


Comparison of all circular, all square, mixed (as in Fig. 1), all diamond and all D-shapes with each other in Table 2 
reveals relatively lower hydrodynamic loads for all four structures observed for the circular shapes at the L/D = 2.00. The 
largest mean drag coefficient here is experienced by the fourth structure in both mixed-shaped and all diamond-shaped 
arrangement. The highest maximum fluctuating drag coefficient of 1.38 is observed for cylinder 2 in the all square-shaped 
arrangement. The largest maximum amplitude of the lift coefficient of 1.63 also belongs to the all square-shaped 
arrangement, but corresponds to cylinder 4. Mixed-shaped arrangement (or basic case of structures with alternate cross- 
sections, as in Fig. 1) at L/D = 2.00 demonstrates a relative consistency in large amplitudes of the lift coefficient for 
cylinders 2, 3, 4, which makes estimations of hydrodynamic loads for this arrangement more important, due to higher 
expected loads. 
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Fig. 2. Time histories of fluctuating drag coefficients for four cylinders of a different cross-sectional shape with L/D = 2.0: 
a — cylinder 1; b — cylinder 2; c — cylinder 3; d— cylinder 4 
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Fig. 3. Time histories of lift coefficients for four cylinders of a different cross-sectional shape with L/D = 2.0: 
a — cylinder 1; b — cylinder 2; c — cylinder 3; d — cylinder 4 


Further observation of signals of the fluctuating drag and lift coefficient in Figures 2 and 3 reveals meaningful instabilities 
in forces experienced by all structures with square shapes, and for structures 2, 3, 4 with diamond and mixed shapes. The 
comparison shows that circular and D-shaped structures would experience lower and more stable fluid loads in this arrangement. 
Table 2 and Figures 2 and 3 confirm that cylinder 1 in the shape of a circle or D-shape provides reducing fluid loads for the 
three downstream cylinders. This substantiates the common interest in further exploration of the effect of the D-shaped structure 
on reducing the drag force in the arrangement with all alternative cross-sectional shapes. 
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This effect is studied in details in the next (second) simulation series, presented for each cylinder in Figures 4—7 in 
terms of the fluid forces and in Figures 8, 9 — in terms of the fluid flow characteristics for the considered computational 
domain. Figure 4 illustrates fluid loads on the upstream structure, where the most unstable signal (at L/D = 2.17) belongs 
to fluctuations of the drag force. Figure 4 c also indicates presence of multiple frequencies in signals of the fluctuating 
drag coefficient, while a single dominating frequency can be identified for the lift force coefficient in Figure 4 d. 
Figures 4 c and 4 d demonstrate that the frequency of both lift and drag forces generally increases with the growing L/D 
ratio, and the maximum frequency is indicated by signals at L/D = 2.67 and 2.83. 
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Fig. 4. Time histories of fluid force coefficients for cylinder 1 in the uniform flow: 
a — fluctuating drag coefficient; b — lift coefficient; c — fluctuating drag coefficient FFT; d — lift coefficient FFT 


Relatively similar complexity of frequencies of the lift force is observed in Figure 6 d for cylinder 3, where one to 
two dominant frequencies could be clearly identified. At the same time, more than two frequencies are observed in 
Figure 6 c for each signal of the fluctuating drag force. The pattern of growth in the overall dominant frequency with the 
increasing L/D is still recognisable for cylinder 3, similar to cylinder 1. Combination of frequencies is even more complex 
for cylinders 2 and 4, as comes from Figures 5 c—d and 7 c-d, this is not possible to indicate clear dependences in the 
frequencies of the fluctuating drag force. Some resemblance of the found growth trend could be still observed in Figure 5 d 
for the lift force coefficient for cylinder 2. This provides evidence for the generally unstable nature of hydrodynamic 
loads acting on the square cross-section shown in Figure 5. 

Impact of the L/D ratio on the mean drag coefficient for cylinder 1, based on Table 3, is partially confirmed: apart 
from a couple of deviations, the mean drag force increases with the growth of L/D. Data for cylinders 2 and 4 do not 
indicate a specific pattern, as the values fluctuate back and forth within 10% from the initial mean drag coefficient at the 
smallest L/D. Cylinder 3, on the contrary, indicates a stable pattern of the reduced mean drag coefficient with the increased 
L/D, so that the reducing ability of the upstream D-shaped structure is evident, but for the central cylinder only. The 
largest mean drag coefficient of 1.02 is linked to cylinder 4 at L/D = 2.67. The considered range of L/D allows observing 
an important transition from the strong to minor interference in the wake of the three paired structures. 

The highest fluctuating drag coefficients of 1.54 and 1.36 are linked to cylinders 2 and 4, respectively, both observed 
at L/D = 1.83. The feature of the maximum amplitude of the fluctuating drag coefficient, indicated in Table 2, is in absence 
of a specific dependence from the L/D, the values rapidly change from near zero to relatively high with a small increment 
of change in the ratio. The fluctuating drag coefficient has generally the lowest amplitudes for cylinder 1, average 
amplitudes — for cylinder 3, and the largest amplitudes — for cylinders 2 and 4. 

The largest maximum amplitude of the lift coefficient occurs at L/D = 2.83 for cylinder 2. The lift coefficient 
generally resembles the distribution, similar to the maximum amplitude found for the fluctuating drag coefficient: the lift 
force appears to be the smallest for cylinder 1, relatively intermediate — for cylinder 3, and the highest — for cylinders 
2 and 4, with no specific pattern linked to the Z/D increase and reduction. This allows us to conclude that the ability to 
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reduce hydrodynamic forces by placing the upstream D-shaped structure in front of the array is limited. The force 
reduction is observed mainly for the structure placed in tandem downstream, and the effect is most pronounced for the 
mean drag coefficient, with some reduced effects also seen for the fluctuating drag and lift coefficients. 
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Fig. 5. Time histories of fluid force coefficients for cylinder 2 in the uniform flow: 
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Fig. 6. Time histories of fluid force coefficients for cylinder 3 in the uniform flow: 
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Fig. 7. Time histories of fluid force coefficients for cylinder 4 in the uniform flow: 
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Fig. 8. Contours of the flow characteristics for L/D = 1.67 at 200 seconds: 
a — velocity contour; b — vorticity contour; c — pressure contour; d — streamline 
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Fig. 9. Contours of the flow characteristics for L/D = 2.83 at 200 seconds: 
a — velocity contour; b — vorticity contour; c — pressure contour; d — streamline 


Figures 8—9 show the velocity, vorticity, pressure and streamlines of the flow around cylinders for some selected L/D, 
where both proximity and wake interference among the cylinders are presented for the time step of 200 s. The flow around 
cylinders is complex, and vortex formation patterns are highly affected, as the distance between the cylinders increases. 
The proximity interference is observed for cylinders 2, 3, and 4, alternate single vortices are shed on the downstream side 
of these structures. For the wake interference at L/D = 1.67, free shear layers separate from the upstream cylinder | and 
reattach themselves to the upstream side of cylinder 3, and a vortex street is only formed at the downstream side of 
cylinder 3. At this distance, a broad region of wake is created at the downstream side of cylinders 2, 3, and 4. As the L/D 
increases above 2.00, there are vortices formed at the upstream side of cylinder 3, in the wake of cylinder 1. A vortex 
street is also formed at the downstream side of cylinder 3, with formation of 2S vortices. Figures 8 a, 9 c demonstrate a 
group of minor vortices formed following the diamond-shaped cylinder 4 and a vortex pair formed in a similar to S+P 
vortex cycle past cylinder 2. 

Discussion and Conclusion. The 2D numerical simulations are performed in this work for cylinders with different 
cross-sectional shapes at the Reynolds number of 2.5x10° using the DES approach. The considered cylinders are studied 
in a complex position of an upstream D-shaped structure in front of three paired structures, with the aim to investigate the 
drag-reducing ability of this specific layout, observe the flow complexity, the wake interference from each structure, and 
vortex formation patterns. 

The following conclusions could be made from this study: 

1. The ability to reduce hydrodynamic forces by placing the upstream D-shaped structure is mainly limited to the 
structure placed in tandem downstream, and the effect is most pronounced for the mean drag coefficient. 

2. Overall, the mean drag coefficient of cylinders is observed to be affected by varying L/D, with the main effect on the 
mean drag coefficient of cylinder 1, which grows with increasing L/D, and of cylinder 3, which reduces with increasing L/D. 

3. Competition of frequencies is observed for the fluctuating drag coefficient for all structures and lift coefficient 
signal of cylinders 2 and 4. This competition is due to the joint effects of both the uniform current and wake interference, 
which intensifies at a lower L/D in terms of changes to the resulting vortex street. 

4. Both proximity and wake interference among the cylinders are observed. The flow around cylinders is complex, 
and vortex formation patterns are highly affected as the distance between the cylinders increased with 2S being the major 
vortex type formed and shedding the additional vortices from the square and diamond structures. 
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The study generally contributes to the field of knowledge by advancing our understanding of fluid-structure 
interactions, drag reduction strategies, and vortex dynamics, with potential applications in offshore energy systems. The 
current work contributes to the development of the drag reduction strategies through analyzing the impact of the upstream 
D-shaped structure on downstream cylinders. Understanding how different structural configurations affect drag can 
inform the design of more efficient systems in various engineering applications, such as subsea transportation of fluids. 
By observing flow complexity, wake interference, and vortex formation patterns, this study contributes to the 
understanding of fluid dynamics around complex geometries. This knowledge is crucial for optimizing the performance 
of structures in environments where fluid flow plays a significant role, such as subsea engineering. The results of the 
present research highlight the effect of varying the aspect ratio L/D on drag coefficients to inform engineering designs of 
similar arrangements. This study reveals the intricate vortex dynamics and shedding patterns, particularly concerning the 
proximity and wake interference among cylinders. Understanding these phenomena can aid in predicting and controlling 
flow behavior around complex configurations, leading to more efficient designs and better performance in practical 
applications in offshore systems. 
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